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Urn models have been studied a lot because they are analytically tractable and contain 
rich physical phenomena, e.g., slow relaxation or condensation phenomena. It has been 
revealed that the urn models are related to zero range processes and asymmetric simple 
exclusion processes^ which are stochastic processes for studying nonequilibrium statistical 
physics. Furthermore, the urn models have been used in research fields of complex networks.'^' ^ 

Recently, the analytical treatment for disordered versions of the urn models have been de- 
veloped.^'®'^ Leuzzi and Ritort have analyzed a disordered urn model (disordered backgammon 
model) in the scheme of the grand canonical ensemble. On the other hand, in the scheme of 
the canonical ensemble, the replica method has been used in order to calculate the occupation 
distribution of the preferential urn model. ^ The replica method is a powerful tool for random 
systems, but still has an ambiguity for the mathematical validation. Furthermore, the analysis 
in ref. 9 has been based on replica symmetric solutions, so that it was not clear that the replica 
symmetric solution is adequate for the disordered urn model, although the solutions are in 
good agreement with numerical experiments. 

In this short note, we calculate the free energy of the disordered urn model using the 
law of large numbers. It is revealed that the saddle point equation obtained by the usage 
of the law of large numbers is the same as that obtained by the replica method in ref. 9. 
Hence, we conclude that the replica symmetric solution is adequate for the disordered urn 
model. Furthermore, we point out the mathematical similarity of free energies between the 
urn models and the Random Field Ising Model (RFIM); this similarity gives an evidence that 
the replica symmetric solution of the urn models is exact. 

Firstly, we give a brief explanation of a general urn model (Fig. 1). An Urn model consists 
of many urns and balls. We consider a system of M balls distributed among N urns; the 
number of balls contained in urn i is denoted by rij, and hence Xli^i ~ '^^^ density of 
the system is defined hy p = M/N. There are two types of urn models: the Ehrenfest class and 
the Monkey class. ^ In the Ehrenfest class, balls within an urn are distinguishable. In contrast, 
the Monkey class has indistinguishable balls. These two types of urn models are treated in 
the similar analytical method, so that we discuss here only the case of the Ehrenfest class. 
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Fig. 1. Illustration of the urn model. Each urn has a disorder parameter hi 



The dynamics of the Ehrenfest urn model is briefly summarized as follows: 

(1) Choose a ball at random. 

(2) Select an urn with a transition probability Ui. 

(3) Transfer the chosen ball to the selected urn i. 

(4) Repeat the above procedures until the system reaches an equilibrium state. 

Selecting an arbitrary transition probability n-j, one can construct various urn models suitable 
for own physical problems. 

The total energy of the whole system is defined as 



N 



(1) 



where E(hi,ni) is an energy of urn i, and hi the disorder parameter for urn i. It is assumed 
that the disorder parameter hi is assigned to each urn by using a probability density (pih). 
The partition function of the urn model is written as^ 
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is the Boltzmann weight attached to urn i, and (3 an inverse temperature. The Kronecker delta 
(5(^^^nj,M) stems from the constraint of the total number of balls. Using the Boltzmann 
weight, the transition probability is defined as^ 

Ui = ^(%^. (4) 

p{hi,ni) 

The aim of the present paper is to calculate the free energy in the thermodynamic limit 
by using the canonical partition function Z of eq. (2). In order to calculate the free energy, 
we rewrite the partition function Z as 

TV 
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where we use the integral representation of the Kronecker delta 

2iTT6{m, n) = ^ dz z^' 
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and define 



H{h,z) = ^p(/i,n)z". 
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Hence, the free energy of the disordered urn model in the configurational average is written 

by 
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where the configurational average is defined as 

{A{hi,--- ,hN)){h} 

= J dhi(t){hi) ■ ■ ■ J dhN(t>ihN)A{hi, - ■ ■ ,/iAr)- 

The law of large numbers is available in order to calculate the configurational aver- 
age of eq. (8). Using the law of large numbers, j^J^iLi^i = (^) ~^ °^)) ^^e term 
{^^iln H {hi, z))/N in eq. (8) is replaced by {In H{h, z))h- Hence, it becomes easy to take 
the configurational average. By using the saddle-point method, we get the free energy 



^{F){f,} = -^\nzs + ^{lnH{h,Zs))^ 



and the saddle-point equation 
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We note that the final result is consistent with the analytical results obtained by the replica 
method.^ Hence, we conclude that the replica symmetric solution in ref. 9 is valid for the 
disordered urn model. 

Finally, we will point out the mathematical similarity between the urn model and the 
RFIM. The RFIM consists of N Ising spins interacting through an infinite ranged exchange 
interaction. The Hamiltonian is given by 



J 



(12) 



i<3 



where J is the spin coupling constant, {Si\ are spin variables, and = {hi, ■ ■ ■ ,/iAr} are 
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the external magnetic fields. The free energy of the RFIM is written as^° 

-^(i^RFIM){h} 




27r T , 
dm exp 

oo 



(13) 

where /? is the inverse temperature. For the RFIM, it has been revealed that the replica 
symmetric solution gives an exact result. We can easily see the similarity of the free energy 
of the RFIM in eq. (13) and that of the urn model in eq. (8): the configurational average 
of a logarithm of an integral is taken, and the term depending on the disorder parameters 
is contained in the integral as the argument of the exponential function. We here remark a 
difference between the urn models and the RFIM. The variable of integration in eq. (13) is 
the magnetization, m = Sj^ /N, and that in the disordered urn model is related to 

the macroscopic constraint M = X^i^i or p = (^^^i rii^ /N. The macroscopic observable 
M (or p) in the disordered urn models is a constraint, therefore the observable is known in 
advance. In contrast, the magnetization m in the RFIM is unknown in advance and constructed 
by microscopic variables Si with time. Although there are such differences between the urn 
models and the RFIM, the macroscopic observable constructed by microscopic variables {Si 
or rii) plays a key role in the both analysis. 

In summary, we discussed a disordered version of a general urn model. The free energy 
and the saddle point equation were obtained by the usage of the law of large numbers, and 
it was clarified that the saddle point equation is the same as that of the replica symmetric 
solution. In addition, we pointed out the similarity between the urn models and the RFIM. 
Although we have not given the rigorous mathematical proof, we expect from the similarity 
that the replica symmetric solution describes exact results in the thermodynamic limit. We 
hope that the analogy with the RFIM will lead to better understanding of the disordered urn 
models. 
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